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ABSTRACT 


Propagation loss due to shadowing by seamounts is studied 
utilizing physical models in air. Dickens Seamount in the 
Gulf of Alaska is approximated by three models: A plane wedge, 
a contoured wedge, and a scaled three dimensional model. The 
forward diffraction over each is analyzed for a five octave 
frequency range. A new concept, the far-field "diffraction 
scattering strength” is defined and used to predict frequency- 
dependent diffraction loss at sea. The total propagation loss 
is calculated by adding laboratory model values of upslope 
forward scatter and crest diffraction losses to computer- 
predicted ray refraction losses up to and away from the seamount. 
This predicted loss is then compared to long range ocean propa- 
gation loss measurements for the case in which rays are com- 


pletely blocked by the Seamount. Close agreement is found. 
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eee RODUCTION 


Diffraction of acoustic energy by wedges has been 
theoretically investigated by a number of authors. Of var- 
ticular interest is the work by Biot and Tolstoy /Ref. Ws 
Their solution to a sound impulse radiated from a point source 
and diffracted by a rigid, infinite wedge employs the approach 
of normal modes and the method of normal coordinates. 

Bremhorst /Ref. 2/7 has demonstrated close agreement 
between Biot-Tolstoy theoretical predictions and experimental 
data using a 90° wedge as the diffracting barrier. 

The question now is whether this predictive ability 
can guide the prediction of forward diffraction by more com- 
Dlex, real-world barriers. Specifically, can diffracted sound 
energy be predicted in the shadow zone of a seamount or under- 
water ridge? 

To accomplish this purpose, three barriers were con- 
structed: a simple plane wedge, a long-crested (two-dimension- 
al) wedge contoured to an at-sea sound track, and a three 
dimensional model of a seamount. Acoustic energy was diffrac- 
ted over each in order to determine "diffraction scattering 
strength", a newly defined concept, and "diffraction loss". 
This was in turn compared to at-sea data showing the effects 
of shadowing on sound propagation by the Dickens Seamount in 
the Gulf of Alaska. The results of this study indicate that 
the diffracted energy over complex barriers can be predicted 
utilizing physical laboratory models. 
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A. ANECHOIC CHAMBER AND DATA PROCESSING EQUIPMENT 

The experiments were conducted in the Naval Postgraduate 
School (air) anechoic chamber, providing superb isolation from 
excess background noise, and a very stable medium for sound 
travel. 

Data acquisition and processing were accomplished using 
a digital computer system composed of four major components 
interfaced to provide high speed analog to digital conversion, 
digital processing and data printout. The design was develop- 
ed by the Special Projects section of Naval Air Develovment 
Center in conjunction with Pinkerton Computer Consultants, Inc. 
of Warminster, Pennsylvania. The major components are: 

1. Interdata Model 79 Computer 
This minicomputer is a digital design that is 
FORTRAN and Basic programmable with a64 thousand byte memory. 
Data that have been stored on floppy discs can be read into 
the computer for processing. There is an identical capabil- 
ity utilizing cassette tapes. 
2. Phoenix Analog to Digital Converters, Model ADC 912 
Two ADC 912 analog to digital converters may be 
used separately or simultaneously. Each converter is a very 
high speed, high accuracy device capable of encoding #10 volt 
input signals into twelve binary bits of data with a 


resolution limit of one part in 4,095 at a maximum time of 





2 microseconds per conversion. It measures the input voltage 
against an internal precision reference voltage with an accur- 
acy of + 0.025% of full scale. The ADC 912's will accomodate 
a typical commutating through-put rate of 476,190 channels 

per second, including settling time. The sampling frequency 
is set in a stable frequency synthesizer and sent to the con= 
verter via a sampling circuit. 


3. Texas Instruments Silent Electronic Data Terminal, 
Model 733 


The TI 733 consists of a key board used as a pro- 
gamming input/output control device, a printer, and a trans- 
mit/receiver mechanism going to a peripheral disc drive unit. 
The overall system facilitates rapid, accurate processing of 
any desired type of analog electrical signal and was used 
primarily for frequency domain analyses using standard FFT 


algorithms. 


B. STANDARD EQUIPMENT LIST 
The following scientific equipment was utilized to con- 
duct the experimental work: 


Interface Technology Timing Simulator/Word Generator, Model 
RS-648 


Tektronix Type 545B Oscilloscope with four trace plug-in 


General Radio Type 1163-a Coherent Decade Frequency Synthe- 
sizer 


Hewlett-Packard Electronic Counter, Model 5223L 
Hewlett-Packard DC Power Supply, Model 721A (two) 
E&L Instruments DD-1 Digi Designer 


Hewlett-Packard Dual-Trace Oscilloscope, Model 140A 
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Fluke True RMS Voltmeter, Model 8920A 

Fluke Digital Multimeter, Model 8000A 

Hewlett-Packard Power Amplifier, Model 467A 

Wavetek Arbitrary Waveform Generator, Model 175 

Lambda Regulated Power Supply 

Krohn-Hite Frequency Filter, Model 3550 (two) 

Princeton Applied Research Amplifier, Mddel 113 (two) 

Bruel and Kjaer Type 4145 Condenser Microphone 

Bruel and Kjaer Type 4134 Condenser Microphone 

Bruel and Kjaer Microphone Power Supply, Model 2804 
The equipment diagram for the experiment is shown in 


Figure 3. 
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Figure 1. Laboratory Computer System 
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Pee cAbORY 


Diffraction can be generally described as a vhenomenon 
which occurs whenever a wavefront encounters an obstacle. 
To examine the process theoretically, consider an infinite 
wedge which 1s bounded by rigid plates at §=0 ; 4 = Be and 
intersecting along the z-axis or apex. The region 
is occupied by a homogeneous compressible fluid of sound vel- 
ocity c and density p é 

To determine the normal modes of the system, it is con- 
venient to solve the acoustic wave equation in cylindrical 


coordinates, in terms of the displacement potential : 


22.1 2, 78 _1 Fe 


me 
oR r or +73 3e a Ct ie) 
The harmonic sokutions to this equation are of the form 
tive), +i, +tiwt 
Q=< Kee de (2) 
where 
2 
2 W) ee 
.? ee (3) 


The rigid plates of the wedge establish boundary con- 


BPeLons such that: 


ot = 0 at B= 0, O= Ow (4) 


Consequently, 


C1, ay + os + iwt 


Q = cos ¥, 8 H e or 


digs 





and AT 
(6) 


At r=0, the imaginary vart of the Hankel function be- 
comes infinite which is not physically realizable. Therefore, 
only the real parts of i) the Bessel functions of the 
meest Kind, will be kept. 

At this point, assume an explosive point source at r=Ior 
G = or , and z=0. This will result in solutions to equation 
(5) which are symmetric in the z-axis. Biot and Tolstoy now 
apply normal coordinate methods. For an in depth look at the 
mathematics of the solution technique, see Refs. 1 and 3. 


Leaping directly to the results, the solution form to equa- 


tion (1) becomes: 


Bape = \i cosy 8 Cos, ©, Ly ctz 
Ot a, 2. : " , 0 (7) 


where = [ T, (Kr) T, Ke) Tk Ct-3°)"]kak (8) 
n n 


A transform given in Appendix X, Ref. 3 may be used to reduce 


I, to more physically visualizeable terms. Letting 


4 r 


| 2 2 
= —[(r- 0) oe) (9) 


6 


t = L (r+) + 3°] /2 (10) 


0 


te is the time required for a sound pulse to travel directly 
Onn 

from source to receiver in the general case. (, is the time 

required to travel from source to wedge apex to receiver, ie., 


the diffracted energy path. As a result, equation (8) will 


have three solutions: 
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—_ “ 
Beret <t , t.< t<(, and (,< t. Since this work is confined 
to studying the energy only in the shadow zone of a barrier, 


only the solution involving the diffracted wave, C< tees! Or 





interest. 
For this case 
= \ ae 
_ 7 su Vite 
ae? ore, sith y (11) 
a 7 . 
ne (Y 4 aa, 
ee Oe 
where eas cosh Zane © (12) 


Now combining with equation (7) 





-D, 4 
2 C ( A A 
ee ) Cos, 6 CHV. 8 sind we Ph <t 
Oot Ta, ryesinh Y 7" ye . ‘s ee ee 


Rewriting the trigonometric functions in terms of their ex- 
ponential identities and then collecting the conjugate pairs 
together, equation (13) becomes ee 

mie 1 SB snCa irtos6) 
Dt Away ye sink y I-22 KicoslG,) rs020,) +6 44, 


For the unit impulse source used in Refs. 1 and 3 the dis- 


(14) 


olacement potential d) is related to the acoustic pressure by 


28 
If, on the other hand, one assumes a source that is a delta 
function in time as well as space, the acoustic pressure due 
to the diffracted wave alone, ie., in the shadow region of a 
wedge 1s given by 7 
-T eee ee +6H+ ) 

et SAE sin )(rto%e 

= SS eB eS ss SS 

. vy sim ae ate ~ 2m 
P Amey Yh Y l-2e 7h, (08 (G, \(r40#6,) + € Va, 





(16) 


ey 








7a. 
where CL = (v4 v4 2») 
Y = are Cosh —-—— 
2YVv 


0 


and 


fluid filled region above the wedge 


a 


o = angle between the wedge side and the source 


Co ®D 
ul 


angle between wedge side and the receiver 


distance from source to wedge apvex 


a 
O 
it 


r = distance from wedge apex to receiver 


zZ = distance out from the least time travel path along 
the wedge apex 


Figure 4 illustrates the typical geometry of the problem. 

The diffraction process may be described in the follow- 
ing manner. A point source at a distance ry from the wedge 
apex is transmitting an expanding spherical wavefront. The 
wavefront first arrives at the apex along the least time trav- 
el path, at time t = r/c The wavefront then continues to 
arrive at progressively later t with increasing distance Zz. 
This has the effect of establishing a series of sources along 
the apex which then reradiate into the shadow zone on the re- 
ceiver side of the wedge. The apex of the wedge is essential- 
ly acting as a line source with time shading along the line. 

Equation (16) has been transformed to the frequency 
domain by Medwin / Ref. 47 and used to compute theoretical 


values which were compared to experimental data from a simple 


wedge / Ref. 2/7. 
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Figure 4. Typical geometry on a simple wedge. 


The source is at r = Tye b=86., and z= 0. 


Receiver is at distance r from the wedqe 
apex at angle . The case oe = 0 and § ate 


describes both source and receiver on the 
wedge surface. 
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Vier POCBAN Gol UDY 


To determine whether the magnitude of the diffracted 
energy could be predicted of course required some real-world 
data to which experimental data could be compared. To this 
end, a study of shadowing of sound propagation by a seamount 
conducted by Ebbenson, et al. /Ref. 5/ was ideally suited. 

A series of propagation runs at a frequency of 230 Hz 
were made by Ebbeson et al. at the Dickens Seamount to deter- 
mine the degree of shadowing and its effect on sound prova- 
gation. To be expected, the sound source depth played an 
important part in the shadowing effect due to its relative 
proximity to a sound channel. But particularly for a shallow 
source where most of the sound rays are deeply refracted, 
there was a significant effect on propagation loss when 
compared to predicted values. With source and receiver 
separation appropriately spaced, it was possible to have a 
major (15 dB) shadowing of acoustic energy when the deeply 
refracted rays totally impinged upon the slope of the seamount. 
See Figure 5 for an illustration of the ray paths. It can be 
seen that at ranges of 79 and 119 km maximum shadowing 
should take place since the source is in a position which 
enables the seamount to intercept the deeply refracted rays. 
Minimum shadowing is expected when the source is at 99 km 


range. 


Figure 6 illustrates the measured propagation loss 


along with that predicted by the FACT ray propagation ray 
20 
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model. This is based on an average sound speed profile shown 
in Figure 7. The measured results agree well with prediction 
up until the point at which the source passes over the apex 
of the seamount (60 km). From there on the increase in 
propagation loss due to the shadowing effect of the seamount 
1S apparent. At A and C there is approximately a 15 dB in- 
crease in propagation loss. 

This increase in propagation loss is here postulated to 
be due to diffraction. The objective of this work is the pre- 
diction of the ocean diffraction loss by physical modeling in 


the laboratory. 
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V. EXPERIMENTAL DESIGN AND PROCEDURES 


A. PHYSICAL MODELING 

To physically model the Dickens Seamount several orders 
of approximation were used. Three models were built, each 
of which progressed to a new stage of complexity. The models 
were required to appear as having infinite dimensions when 
compared to the wave length of the sound signal. This was 
done by ensuring that the edges of the model are sufficiently 
far away from the receiver so that a single cycle pulse at 
the lowest frequency could be received without interference 
from edge diffraction. Therefore, with sufficiently large model 
dimensions in conjunction with a pulsed signal technique, the 
models were considered infinite for the frequency band in 
which the experiment was conducted. Furthermore, it was desired 
that the model surfaces approximate a perfectly rigid condition. 

1. The 14° Wedge 

The first model and the simplest was constructed on 

the assumption that the seamount behaves acoustically as if 
it were a simple, geometric wedge. That is, for this model 
it is assumed for simplicity that the peripheral portions of 
the seamount make little contribution to the diffraction 
process and the major effects are along or close to the least 
time travel path. Additionally, along that travel path there 


1s an average plane upslope and plane downslope. 


ZS 





Therefore, a simple plane wedge was constructed to 
correspond to a particular track across the seamount. Re- 
ferring to Figure 8, this is track 6 of the ocean exveriment 
/Ref. 57. The average upslope and downslope measured with 
the horizontal is 14° so this is approvriately named the "14° 
wedge". It was made of two pieces of 1/4" thick aluminum 
each 1.52 m (5 ft.) long by 0.6 m (2 ft.) wide and joined at 
the apex to make a 152° interior angle con = 212° in the no- 
tation of Ref. 1). Aluminum was chosen as the building mater- 
ial due to its rigidity and its reflection coefficient of 
almost unity. 

2. The Contour Wedge 

The next level in modeling complexity was to go 
a step beyond the concept of an average slope and to account 
for a contour. Therefore, the wedge was given a smoothly 
varying contour approximation to that of the actual seamount 
along track 6. Acoustically, the wedge still appeared in- 
Finite but the wavefront had to interact with a variable 
sloped surface rather than a flat plane on each side of the 
wedge apex. The contour did not vary with respect to the 
Z-axis. 

To construct this "contour wedge", as it was 
called, the technique resembled that used in aircraft wings. 
Four 1" x 12" pieces of lumber were cut conforming to the 
track 6 contour. These in turn were fastened by their flat 
Sides to a 3/4" plywood base. Then .040" thick sheet 


aluminum was laid over the 1" x 12" formers and fastened with 
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TRACK 6 CONTOUR 


Contour mav of the modeled area of Dickens seamount. 
Paths 1 through 6 represent least time travel vaths 
Ceeeeomenmcze km coacing utilized in the laboratory 
acoustical survey. 
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countersunk screws making the surface of the contour wedge. 
Aluminum was again used for its reflection coefficient and 
workability. The underside of the aluminum surface had to 
be coated with ashpalt roofing compound to add rigidity and 
dampen out unwanted vibration in the presence of sound. The 
contour wedge was 1.83 m (6 ft.) wide by 1.22 m (4 ft.) long 
(in the z-axis). 

3. The Three Dimensional Model 

The final stage of complexity was a three dimen- 
Sional model of the seamount. The Defense Research Establish- 
ment Pacific (Victoria, B.C., Canada) provided superbly 
detailed bathymetric data from which the 3D model of Dickens 
was constructed. 

The at-sea experiment was done using a frequency 
of 230 Hz. To duplicate the experiment in air required either 
higher frequencies than the A/D converter could handle or a 
larger model size than could be accomodated by the entrance 
door to the anechoic chamber. As a compromise, the 3D model 
base was made approximately 2.13 m (7 ft.) square. The 
scale was 1/2" = 100 m. Because of the change of medium 
this is equivalent to 1:7874. When irradiated with sound 
of frequency 100 kHz in air, the model represented a fre- 
quency of 55 Hz. in water at Dickens. 

The 3D model was constructed by laminating 1/8" 
layers of particle board onto a 1/2" particle board base. 
Each layer was cut from a properly scaled blow-up drawing 


of the bathymetric data with constant depth contour lines 


28 





p 


a4 


ve 


Ou, 


6 


DANHT 





Zo 





SADE MMAR 


iY 


of: 


Gabe: 
4, 


py 
Ze 
LO 
‘4, 


t 


A 
COPE Ly 
MODES, 
Woes 





30 


The Contour Wedge 


Fuguiee U0: 





T=pol 


q 


ele eae 


qd 


€ ouL “TT 





PAA 


NTRNAAAA BA Sa 


FLLHON 


ha ae a ated 








to a depth of 2500 m marked in. This technique resulted 

in a depth resolution on the 3D model equivalent to 25 m 

at sea. Once the layering was complete, the entire model 
surface was coated with a thin layer of plaster to fill in 
the steps between each layer of particle board and to give 
the surface a smooth appearance. The combination of pvar- 
ticle board and plaster was of sufficient density to result 
in an acoustically rigid reflector. The builder was Model- 


makers, Inc. of San Francisco. 


ge SOURCE/RECEIVER SELECTION 

To fulfill the conditions of the Biot-Tolstoy Wedge 
Theory which guided this work, it is desirable to have an 
approximation of a point source to transmit spherical waves. 
This requires that the source must be small compared to the 
wavelength of the transmitted signal. This condition is met 
if ka<< 1 where k is the wave number and a is the radius 
of the sound source. Additionally it is required that the 
radial distance from the source to the apex of the barrier 
be much greater than X. 

With the above in mind, the selection of a sound source 
becomes one of compromise. To achieve the closest approxi- 
mation to a point source requires the smallest possible 
radiator; however, this also limits the maximum acoustic 
pressure that can be generated. 

Bremhorst was confronted with the same problem of source 


selection and ran extensive studies to determine the opvtimum 





radiator subject to the above constraints. He determined 
that the best solution was to utilize the reciprocal trans- 
ducer properties of small condenser microphones and employ 
them as sound sources. 

Likewise for this work, Bruel and Kjaer (B&K&) condenser 
microphones were put to the task of acting as the sound source. 
These microphones, although being piston radiators, do have 
the fundamentally important property of radiating a spherical 
diverging wave but have a directivity term involved. This 
directivity factor is rendered insignificant by aiming the 
microphone at the barrier apex and ensuring the maximum re- 
sponse axis is along or very close to the least time travel 
path between source and receiver. 

Laboratory testing was performed using the 1", 1/2", 
and 1/4" B&K microphones. It was found that only the 1" B&K, 
Type 4145, had a sufficient acoustic pressure output at 190 kHz 
to be used. This requirement for the maximum acoustic pressure 
attainable was particularly evident when transmitting sound 
over the 3D model because of the undulations along the sound 
path which caused some secondary scatter adding to the trans- 
mission loss. See Figure 12 for an illustration of the theo- 
retical acoustic output of the various B&K microphones. 

The selection of a receiver was again a matter of com- 
promise. In this case the trade off was between trying to 
use the smallest microphone possible to ensure good high 
frequency response but at the same time preserving a suffi- 


ciently high sensitivity level. It was determined that the 
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1/2" B&K, Type 4134, was the best choice, with selection being 
Drimarily based on a useable sensitivity level. Experiments 
were also run to determine if a 4 mm I.D. probe could also be 
used with the microphone to improve high frequency performance. 
Based on typical pulse lengths of 400 microseconds, a 7 cm 
probe was required to prevent interference from reflected 
energy of the probe tip. Unfortunately, the probe resulted 
in too great a loss in acoustic pressure reaching the receiving 
microphone and therefore could not be used. 

In the case of both the source and receiver, they were 
of sufficiently small size to still meet the requirement of 
approximately a point source, at least for the lower fre- 


quencies utilized in the experiment. 


. SIGNAL PROCESSING 
Pe SOUnCe joomla L 

The transmitted acoustical signal waveform must 
consist of many frequencies so that the analysis of received 
data can cover a broad spectrum in the frequency domain. This 
was achieved by using a triangular waveform. It has a sharp 
positive going pulse which acts as an impulse and, when 
transformed to the frequency domain, provides frequencies at 
every harmonic of the fundamental. The harmonic amplitude 
Spectrum is at 6 dB per octave. 

To generate this waveform the Wavetek Model 175 
Arbitrary Waveform Generator was utilized. With the ability 


to program the desired waveform plus a rise-time of 509 ns, 


BS 





it is an ideal frequency generator for the experiment. The 
short rise time makes a good approximation of the infinite 
slope of an impulse. 

The output signal of the waveform generator was 
pulsed to provide the ability to selectively sample only the 
diffracted signal and to eliminate any interference caused by 
reflections from surrounding structures. Data were collected 
using one cycle of the triangular waveform. Triggering of 
the waveform generator to provide each sent pulse was ac- 
complished by the Interface Technology Model RS-648 timing 
Simulator. 

2. Received Signal Processing 

After the diffracted acoustic siqnal was received 
by the 1/2" B&K microphone, it was amplified, band pass fil- 
tered to eliminate low frequency noise and prevent aliasing, 
and then amplified again for a total gain of 46 dB. The 
amplifiers used were PAR model 113 preamps due to their very 
low self-noise characteristics. This signal was then sent 
to the Phoenix A/D converter for digitizing. 

To separate the diffracted signal, a samoling 
Window was used. This was accomplished by the Interface 
Technology timing simulator which triggered open a sampling 
circuit to coincide with the beginning of the diffracted 
Signal. This sampling window could then be left open for 
any predetermined time, then shut. This technique gates out 


any unwanted signal with 100 nanosecond accuracy. 
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The sampling curcuit was built primarily with 
IC devices, associated power supplies, and an oscillator. 
A schematic of the circuit is illustrated in Figure 13. A 
general Radio Model 1163-A Decade Frequency Synthesizer was 
used to deliver the sampling frequency. Since it has the 
Stability characteristics of a crystal oscillator, it pro- 
vided a constant frequency to the A/D converter thereby in- 
creasing the accuracy and repeatibility of the data analysis. 
The sampling circuit works in the following manner. The 
frequency synthesizer puts out a constant frequency sine wave. 
The sine wave is fed to a LM710CN voltage comparator which 
Changes it to a square wave required by the A/D converter. 
Since the frequency synthesizer was not triggerable, a 
. 74123N retriggerable monostable multi-vibrator was used to gate 
open and closed the sampling window. The timing simulator 
provided the trigger signal, coinciding with the beginning 
of the diffracted signal, to bring the "one shot" to a high, 
or on state. By then adjusting a variable resistor in series 
with a capacitor, the high state could be extended to the 
Sesared length before dropping back to low. This procedure 
established the sampling window. Simultaneously, the high 
state of the "one shot" and the square wave output of the 
voltage comparator were fed to a DM5411N AND gate. The output 
of the AND gate was a train of rectangular pulses, only for 
the duration of the sampling window, and sent to the A/D 


converter. These pulses in turn triggered the A/D converter 
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to sample the analog signal. Each rectangular pulse corre- 
sponded to one-half the period of the sampling frequency 
set on the frequency synthesizer. 

The number of samples gathered by the A/D converter 
in each sampling window was set to be a power of two, i.e., 
64, 128, 256, etc. This was required by the FOURONE FFT 
algorithm which was being used to do the frequency domain 
analysis. The sampling frequency was set at an integral 
multiple of the frequency of the sent signal. Due to the 
excellent frequency stability characteristics of the fre- 
quency synthesizer and the Wavetek 175, this relationship 
could be precisely controlled thereby minimizing any trun- 
cation of the sampled waveform and eliminating sidelobes. 

Re wusenie relma rove Cie iqudlity Of the information 
being digitized and submitted for Fourier analysis, an averag- 
ing routine, performed by software, was employed. The data 
from as many as 9,999 samoling windows could be averaged 
prior to performing the FFT; however, only 1,000 windows 
were normally used. This improved the signal-to-noise ratio 
by 10 log N where N is the number of windows averaged. If 
N = 1,000, this means a 30 dB improvement in signal-to-noise 
aoa ©. 

The experiment was originally begun using a 2.5 kHZ 
fundamental triangular waveform as the transmitted signal. On 
the receiving end, 1,000 windows, each of 400 microseconds 


duration, and each containing 128 samples were time averaged 


So 





prior to submission to the FFT algorithm. This resulted in 
Satisfactory data on the diffracted sound behavior up to 
approximately 70 kHz. However, as indicated earlier, this 
behavior must be extrapolated to 430 kHz. To get more 
energy into the higher frequencies, it was decided to run the 
sound source at near its resonance frequency which is approxi- 
mately l0kHz. Consequently, the experiment was continued 
using a 10 kHz fundamental triangular waveform which extend- 


ed the satisfactory results to 100 kHz. 
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VI. EXPERIMENTAL RESULTS AND ANALYSIS 
A. DIFFRACTION SCATTERING STRENGTH (DSS) AND DIFFRACTION 

LOSS (DL) 

In dealing with underwater acoustic behavior it has 
always been convenient to specify the scattering character- 
istics of objects in decibels. For example, in the sonar 
literature "target strength" (TS) is used to give a decibel 
measure of the amount of backscatter assuming spherical 
radiation from the source and spherical reradiation from the 


target. In algebraic terms 


Ens R* 
Pp RR 
oO oe | 


TS = 20 log (17) 


where 


rg 
i 


backscattered acoustic pressure 


J 
it 


reference pressure 
R, = a = reference distances 


For the case of diffraction it is desirable to apply a 
Similar universal concept. If one considers the reradiated 
acoustic energy from a wedge as scatter, then it is appropri- 
ate to define the decibel measure of that scatter as "dif- 
fraction scattering strength" (DSS). 
Clearly, equation (17) must be modified by changing Pic 
to a term which describes the pressure due to the diffracted 


wave alone, i.ée., P(0, 9, On ' har) . With respect to the 


4l 





work being reported here, it is also necessary to be able to 
Scale the DSS from water to experiments in air. To accomplish 


this Medwin /Ref. 4/7 defines the DSS to be: 


TSS = 20 log] Plen OG for) 2 \< 


6 
where r=rfr 


R_ = 1M reference distance 

‘= Y% 
Equation (18) is dimensionless and accounts for spherical 
divergence from a second source, followed by cylindrical 
divergence after diffraction over the wedge apex. Furthermore, 
for a wedge, the equation gives the same value for air or 
water. It has been shown /Ref. 4/7 that the DSS becomes only 
a function of Cuno and &,, provided the ranges are equal and 
sufficiently large, i.e., r = a A : 

The decibel measure of "diffraction loss" may be de- 

fined as: | 


DL =-2o0 | os —20 oo] | 


(19) 


where Pp = rererence psessure and ta is direct distance to receiver. 


From equation (18) 


P Oy 0, Ow Yow") 





Goegires 
MS — 20 eg he JE ee 


Therefore; the diffraction loss becomes 


20 loq 
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oni fq 
Peas 2 O log BE - 20 YR, (21) 


For the results which follow, equations (18) and (21), 
Giving diffraction scattering strength and diffraction loss 


respectively, are applied to the experimental data. 


Be COMPARATIVE DIFFRACTION BY THE MODELS 

The diffracted pressure fields in the shadow regions of 
the three models were studied first to quantify their dif- 
ferences. The 14° wedge provided baseline data with which 
data from tne more comolex models could be commared. Also, 
it provided a quick experimental verification of the Biot- 
Tolstoy theory transformed to frequency space for a shallow 
sloped wedge. 


With the source positioned at 8, = 0 and ae 25), So and 
45 cm, the receiver distance r was keot equal to r5 and O= 0,)- 


On the 3D model, the source/receiver position was along 
track 6 /Ref. bye The decibel measure of the diffracted acous- 
tic pressure was used to compute the diffraction scattering 
strength. These sepvarate results are illustrated graphically 
in Figure 14 through 22. The theoretical DSS, determined 
from the Biot-Tolstoy theory, is also plotted along with the 
experimental DSS of the 14° wedge and shows very good agreement. 

Figure 23 gives a good synopsis of the results. Note 


that the DSS for the 14° wedge remains at a relatively constant 
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-1 dB throughout the range of frequencies studied. This may 
be attributed to the fact that only the single, sharp crest 
of the simple wedge is doing all the diffracting for all fre- 
quencies. The £% behavior of the Biot-Tolstoy theory is 
effectively absorbed by the definition of DSS. There are no 
other variations along the sound travel path for the wave- 
fronts to interact with. 

Sound propagating over the contour wedge, however, is 
confronted both with a rounded rather than sharp crest and 
with a gently undulating surface rather than a plane surface. 
The variations introduce secondary scatter thereby lowering 
the DSS. Nevertheless an asymptQtic value is reached for 
Vs, > 50. 

The 3D model, with a much more irregular surface, offers 
additional secondary scattering effects. Again, the secondary 
scatter is frequency dependent; as the wavelength becomes 
comparable to or smaller than the various dips and bumps in 
the contour, secondary scatter takes place. At n = 1o 
secondary scatter is at a minimum and the DSS of all models 
approach the same value. The frequency dependent behavior 
is clearly evident for 10 ¢ Za SP ao we An Wa OER ere 
is again an asymptotic, far-field, value, although it is lower 
than for either of the simpler models. 

These results clearly show that the experimental study 
of diffraction by complex barriers, and particularly any 


attempt to duplicate or predict real-world data, requires a 
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three dimensional (3D) model. The remainder of this work 
will be based on data taken from the 3D model of Dickens 


Seamount. 


SA PeOUST ICAL SURVEY OF THE 3D MODEL 
To further understand the variability in diffracted 
acoustic energy by the 3D model, an acoustical survey was 
verformed. 
[eeoehSleLvity To Contour 

The sensitivity to contour along the sound travel 
path (track 6) was investigated by varying 6 and r to 
successively greater symmetrical distances from the crest. 
The angles of the source and receiver were kept at 9, = 0, 
O=6 . Figure 24 illustrates the results. As Io and r are 
moved, the falloff of the DSS is different as more of tne 
model contour is acoustically exposed thereby activating more 
secondary scatter points. At A > Ope ome BCC Meir Ves inchs 
into the same asymptote, DSS = -12 dB. 

These data also show that the DSS is apvoroximately 
constant at a given A regardless of whether r or A is 
varied. 

Looking at Vee 50 there is evidence that the 
diffracted energy reaching the microphone is a summation from 
the various radiation points with some degree of coherence. 
In the case studies, an average curve coulc be given such 
that the secondary scatterers result in a variation of + 2 dB 


yr 
around the average for Vee 50. 
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2. Sensitivity To Position Of Ridge Crossing 


TieeselsterwaicycOethe DOsLtion of ridge crossing 
was investigated by transmitting sound over the model at 
various paths roughly perpendicular to the ridge line. Figure 
8 depicts the path locations on the seamount. The source/ 
receiver distance was kept constant at 45 cm (3.54 km) and 
2) = 0, O= Oi 1.e€., source and receiver on the surface. 
Figures 25 through 30 illustrate the sevarate results. 

Figure 31 is a composite of the experimental curves. Figure 

32 gives a comparison of contours along each path. Note that 
there is no great change in ti and consequently there is little 
difference in DSS at A = 10 for the various paths. At high- 
er frequencies, however, the distinct contours each have their 
own secondary scatter characteristics and the different DSS 
curves reflect this. 

Figure 25 warrants some additional comment. As can 
be seen, the DSS curve is atypical compared with other data. 
The 3D model undergoes some radical changes in topography close 
to path 1 (See Figure 8). To one side of the path there is a 
rapid falloff in elevation. It is believed that this extreme 
change in contour is reseosi bie for interfering secondary 
scatter, accounting for the variability in DSS when it should 
be asymptotic. 

It also should be noted that for the duration of 
the data taking for this section, there was a noticeable in- 
crease in electronic interference within the building resulting, 


in particular, in great variability of the 88 kHz data point. 
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3. Sensitivity To Sample Window Duration 


The DSS also has a dependence on how much of the 
model has been acoustically exposed during the opening of the 
sampling window. To investigate this sensitivity, three 
sampling window durations were used: 102, 204, and 410 as. 
These corresponded to receiving energy being reradiated by 
the ridge of the seamount, at distances 11.2 cm (882 m), 
iemeeem (1.76 km), and 23.5 cm (3.53 km, respectively, measur- 
ed on both sides of the least time travel path. In theory, 
as the sampling window is extended, the DSS should increase to 
a maximum value and then remain constant. This is due to 
relatively constanc& reradiation of energy from the wedge crest. 
out to approximately 20 cm (1.6 km) from the least time travel 
path or z = 0 /Ref. 2/7. Beyond that, the falloff of reradiated 
energy is rapid and the contribution to the diffracted pnulse is 
negligible. 

Figure 33 illustrates the results of this investi- 
gation of the 3D model. As can be seen, when the sSamoling win- 
dow is extended from 192 to 204 US , the DSS at the higher 
frequencies decreases. It is surmised that partial destruc- 
tive interference is taking place at the greater window open- 
ings. However, the difference in DSS using the 204 and 410 ys 
sampling window is only about 1 dB, suggesting that at the 
larger openings the interference has reached its maximum and 
that contributions from the greater ridge distance as predict- 


ed are no longer significant. 
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4. penst El yEey Bor AZ ma eh 


To investigate the sensitivity to the angle of 
ridge crossing, or azimuth, the source was set at ie aa 7em, 
0. = Q, and z5 = O with the least time travel vath in the center 
of the model. The receiver was positioned at r = 45 cm, 0=6. ’ 
@ieeaz was Varied from z = 0 to z = +t 25 cm, + 50 cm. These dis- 
tances corresponded to + 2 km and + 4 km at Dickens. Figure 34 
illustrates the results. At angles off from the least time path 


the diffraction loss is greater. 


D. COMPARISON OF LABORATORY PREDICTIONS WITH AT=$SEA DATA 

The total propagation loss measured at Dickens for the 
case where all the rays from the sound source are blocked by 
the seamount is approximately +100 dB re 1M (Case A, Fig. 5). 
It is postulated that this propagation loss is composed of four 
distinct components: (1) the propagation loss (spreading and 
refraction loss) from the source to the side of the seamount, 
(2) a forward scattering loss when the wavefronts imoinge upon 
the seamount surface and propagate toward the crest, (3) dif- 
fraction loss going over the seamount, and (4) the propagation 
loss from the crest of the seamount to the receiver. If these 
separate propagation losses can be determined analytically or 
experimentally then, if the diffraction description is correct, 
their summation should approach what was actually measured at 
sea. 

The prediction therefore combines two approaches: fray 


theory to describe the frequency independent refraction from 
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source to seamount to receiver, and physical modeling based 

on wave theory to describe the frequency dependent provaga- 
tion over the seamount. For this work the propagation losses 
due to refraction and divergence were determined by computer 
at the Naval Postgraduate School utilizing FACT models which 
were parts of the Integrated Command ASW Prediction System 
(ICAPS). The forward scattering and diffraction losses at 

the seamount were determined experimentally using the 3D model. 

1. Measurement Of Scattering Loss 

As the wavefronts emanating from the source con- 
tact the surface, forward scatter takes place which directs 
energy upslope to the crest and parallel to the surface. 
Theories are available for the prediction of forward scatter 
from a statistically rough plane surface. However, no in- 
formation about the surface roughness at Dickens is available. 
Furthermore, the topographical roughness of the upslope of 
Dickens along track 6 is what 1s specifically needed. There- 
fore, an experimental measurement of the topographical forward 
scatter was performed. The objective was to compare scattered 
acoustic energy near the crest with incident energy upon the 
surface. 

The low Q sound source used for the particular 
exveriment was a 8 by 9 cm rectangular solid dielectric trans- 
ducer. The measured radiation patterns for 50 and 100 kHz 
are shown in Figure 35. The receiver was a 1/2" B&K micro- 
nhone because its size was small enough to collect scattered 


Peaatation at the surface. 
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Referring to Figure 5 it can be seen that the 
middle of the cone of incoming rays from the source corre- 
sponds closely with the midpoint of the upslove. Consequent- 
ly, this was chosen as the aim point for the maximum response 
axis of the source. The receiver was placed at a position 
On the upslope corresponding to the intercept point of the 
uppermost rays. The geometry of the experiment is illustrat- 
ed in Figure 36. 

The transient signal utilized was a single cycle 
of a 50 kHz and then a 100 kHz sine wave. The results are 
given in Figures 37 and 38. The first oscillation is the 
direct path signal; this is followed by 3 or 4 cycles of 
receiver ringing. The arrival of the major portion of the 
scattered radiation is indicated by the large amplitude 
oscillations at approximately 120 Ss following the direct 
oath arrival. 

To estimate the decibel measure of the scatter- 
ing loss the ratio of the voltage amplitude of the two 
Signals is used. For both the 50 and 100 kHz cases the 
scattered signal is approximately double the direct path 
Signal. Corrections are needed for the relative distance 


(B) and the off axis directionality (C). Therefore 


SCATTERING 


Abas GC + § db (22) 
LOSS 


dB 
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Figure 35. Beam pattern for dielectric transducer 
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Figure 36. Geometry of the forward scatter measurement 
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where 





VOIRECT 
A = 20 LOG 
VSCAT TERED 
R OIRECT 
B = 20 LOG 
SeAmteReED 
C = dB DOWN AT 5.5° 


Eeg., scattering loss at 50 kHz = -4.8 -2.8 -1.5 +6.0 = 

-3.1 dB. A correction of +6 dB is required since there is a 
pressure doubling at the rigid boundary of the model surface 
which is part of the scattered energy. This results in a for- 
ward scattering gain of approximately +3.1 aB at 50 KHz and 

+ 6.3 dB at 100 kHz. 

It should be recognized that using a pulsed signal 
results in a "snapshot" of the scattered signal. To get an 
accurate measure of the scattered energy using CW at sea, a 
CW signal should be used in the Jaboratory. This would in 
effect integrate the reverberation over a suitable time. It 
is presumed that the peak forward scatter found using the 
pulsed signal is close to the true scatter. A comouter study 
is presently underway to resolve this problem, but the results 
were not available at the time of this writing. 

The degree of forward scatter here is a function of 


what could be called the mesoscale roughness. The modeled 
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FIiguLe 37. Oserloscope disolay Of direct 
path and forward scattered signal; lv/cn, 
100Ms/cm, Ff = 50 kHz. 





Figure 38. Oscilloscope display of direct 
path and forward scattered signal; lv/cm, 
100 yg /em, £ = 100 kHz. 
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surface reproduced the bathymetric readings to increments 
of 25M of deoth in the ocean. This is equivalent to 1/8" 
On the model and is approximately the wavelength in air of 
a 100 kHz signal. It is assumed that the scattered energy 
is not particularly sensitive to the unknown fine, local 
roughness of order less than one wavelength at the intercent 
point on the seamount; but rather it is most affected by the 
mesoscale roughness along the entire travel path along the 
Ssilepe. 
2. Measurement Of Diffraction Loss 

It 1s the seamount surface-scattered energy, 
traveling over the upslope surface, that diffracts over the 
crest. Therefore, the at-sea situation was experimentally 
modeled with the source on the surface of the 3D model. 

Receiver orientation was based on those rays 
that arrive at the receiver after diffraction. Since the 
DSS is strongly dependent on 6, this is a critical decision. 
This orientation was determined by making the assumption that 
the crest of the seamount acts as a reradiation line for the 
diffracted energy. By designating the crest as a sound 
source, a ray trace program was run by Fleet Numerical Weather 
Central in Monterey to ascertain what rays arrive at the 
receiver. From Figure 39 it can be seen that only two rays 
intercept the receiver at the 329M depth at 60 km range: 
the +2° and -6° rays measured with respect to the horizontal 


at the crest. 
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Figure 39. Ray tracings with seamount crest 
designated at the source. 
Receiver denth = 329 m. 
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The diffraction loss experiment was set up with 


a x = 45 cm (3.45 km at sea), which is close to the mid- 


oO 


point of the upslope; 6 = 0, and 6 = +2°, -6°, and -14° 


oO 
with respect to the horizontal. Figure 40 illustrates the 
resultant DSS. The DSS for 0 = +2° is asymptotic at 
-3.0+0.5 dB and for 0 = -6° the asvnpeeteurse—o,020,5 ds. 

As can be seen the DSS is constant for “A PaeeAole 
At the upper ray intercept point to crest range, ry at sea, 
we have A = 543. Therefore, utilizing both DSS asymptotes 
Memeatacion 21, the diffraction loss is +24+0.5 and +2710.5 dB. 
for the +2° and -6° rays respectively, r._ = r = 3.54 km at 


Oo 
2390 Hz. 


3. Computation Of The Refraction Losses 

The first propagation loss computed using the 
FACT model was the refraction loss between source and seamount 
surface at the point where ry = 3.54 km the individual rays 
intercepting the seamount were run through the program. The 
average propagation loss was found to be +76 dB re 1M with a 
+ 2 dB variation depending on the ray angle. 

The computation of the additional refraction loss 
from seamount to receiver was complicated by the fact that 
the +2° and -6° rays from the crest do not originate at the 
source reference position for propagation loss. To solve 
this problem, oe sound source identified in the computer 
program was set at the correct range (19 km from the crest) 
but displaced in depth to allow the generation of grazing 


rays of +2° ana -6° SUETOLNGeErEom ene Crest. The difference 
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in propagation loss to a range deyond the crest equal to r 
(3.54 km) and to the receiver is then the required additional 
propagation loss necessary to account for the total vropna- 
gation loss on both sides of the seamount. By this method 
the additional propagation loss was determined to be +6 dB 
for the +2° ray and +12 for the =G2 ray. 
4. Total Propagation Loss 
The contributing propagation losses are given 


in the following table: 


oe ray -6° ray 

Propagation loss from source +76.0 + 2 dB 6 Out 2. ais 
to seamount upslope at 

r= 3.54 km (for the range 

of incident rays) 
Rough surface scattering -6.3 dB =62.5 a5 

loss at upslope 
Diffraction loss over PaO tee oc te at 0. 9 aS 
crest at oe 3.54 km 
Additional propagation loss +6.0 dB tears 0 aS 

from r = 3.54 km beyond 

crest to receiver 
Total Propagation Loss Po seatoeneoudaet LOO.3 £ 2.5 dB 


The contribution of the two arriving rays at the 
receiver may be summed, on the assumption that the energy 1s 
incoherent, to vield the total predicted propagation loss of 
memes 2.5 GB re IM. For comparison, from Figure 6 the at-sea 
propagation loss in the fully blocked case is +98 + 5 dB re 


r 
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Vinge CONCLUSIONS 


Experimental investigation of the nature of diffraction 
over complex barriers can be accomvolished in the laboratory 
using physical models. 

The new concept of dimensionless diffraction scattering 
strength has oveen shown to be effective in determining range 
and frequency dependent diffraction loss from a scale model. 

The laboratory prediction of long range propagation 
loss was found to be within 1 dB of the propagation loss 
measured by Ebbeson, et al., at Dickens Seamount. With proper 
mesoscale surface detailing of three dimensional models, it 
has been demonstrated that accurate predictions of at-sea 


shadowing loss can be obtained. 
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